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Abstract 

Directed spiral percolation (DSP) is a new percolation model with crossed external bias fields. 
Since percolation is a model of disorder, the effect of external bias fields on the properties of 
disordered systems can be studied numerically using DSP. In DSP, the bias fields are an in-plane 
directional field (E) and a field of rotational nature (B) applied perpendicular to the plane of the 
lattice. The critical properties of DSP clusters are studied here varying the direction of E field 
and intensities of both E and B fields in 2 dimensions. The system shows interesting and unusual 
critical behaviour at the percolation threshold. Not only the universality class of DSP model is 
found to belong in a new universality class than that of other percolation models but also the 
universality class remains invariant under the variation of E field direction. Varying the intensities 
of the E and B fields, a crossover from DSP to other percolation models has been studied. A phase 
diagram of the percolation models is obtained as a function of intensities of the bias fields E and 
B. 

Keywords: Disordered systems, percolation, external bias fields, anisotropy, critical exponents 
and universality. 

1 INTRODUCTION 

Recently, a new site percolation model, directed spiral percolation (DSP)^H], was constructed 
imposing both the directional and rotational constraints simultaneously on the ordinary percolation 
(OP) 3 model. A new type of anisotropic percolation cluster, different from directed percolation 
(DP)Hj and spiral percolation (SP)[3], was found in the DSP model. The model showed unusual 
and interesting critical behaviour at the percolation threshold for the clusters as well as their hulls. 
It was found that the DSP model belongs to a new universality class than that of other percolation 
models. 

There are certain physical properties of a system needed to be studied in the presence of external 
biases like electric or magnetic fields. The electric field gives rise to a directional constraint and the 
magnetic field gives rise to a rotational constraint on the motion of a classical charged particle in a 
plane when the magnetic field is applied perpendicular to the plane of motion. Since percolation is 
a model of disorder, Hall effect in disordered systems can be studied using DSP model. Motion of 
charged particles under crossed electric and magnetic fields in disordered systems is an important 
field of study in recent time UJ. 

In this paper, the DSP model is studied on the square and triangular lattices in 2-dimensions 
{2D) varying the field intensities and the direction of the fields. As the direction of a field changes, 
the number of components of the field on a lattice also changes. The universality class of DSP 
model has been verified for different number of field components. It is found that the universality 
class of DSP is independent of the field directions or the number of components of the directed 
field. As the intensities of the fields change, the DSP model is expected to change to OP, DP or 
SP models at the appropriate field intensities. Varying the external field intensities, crossover from 
DSP to other percolation models has been studied. A complete phase diagram is obtained for the 
percolation models studying the DSP model under variable field intensities. 

Below, the DSP model will be developed under variable field directions and variable field in- 
tensities. The effect of change of field direction will be discussed first. Intensity effect will be 
considered later. 



2 DSP Model under Variable Field Direction 



The DSP model for different field orientations is constructed on the square and triangular lattices 
in 2D. The particles are considered to be classical charged particles to form percolation clusters 
under external bias fields. Two external bias fields perpendicular to each other are present in the 
model. An in-plane electric field E gives rise to a directional constraint and a magnetic field B 
perpendicular to the plane of the lattice, directed into the plane, gives rise to clockwise rotational 
constraint to the occupation of a lattice site by a particle. A single cluster growth Monte Carlo (MC) 
algorithm to generate DSP clusters under different field orientations will be described here. In this 
algorithm, the central site of the lattice is occupied with unit probability. All the nearest neighbors 
of the central site can be occupied with equal probability p in the first time step. As soon as a 
site is occupied, the direction from which it was occupied is assigned to it. Once the directional 
E field is fixed, the empty nearest neighbours in that fixed direction in space are accessible to 
occupation due to the directional constraint. In contrary, for a given B field direction, the empty 
nearest neighbours accessible to occupation are in the forward direction or in a rotational direction, 
clockwise here, with respect to the direction of occupation of the present site. Since the forward and 
rotational directions depend on the direction of occupation of the present site, then the rotational 
constraint is not fixed in space. First, a list of empty nearest neighbours eligible for occupation in 
the next MC time steps is prepared. Selection of eligible sites on the square lattice is illustrated 
in Figs^a) and^6) and on the triangular lattice it is illustrated in Figs^c) and^a!). On the 
square lattice, the E field is in the horizontal direction in Fig^a) and it is along the upper left to 
the lower right diagonal of the lattice in Fig^fc). On the triangular lattice, the E field is horizontal 
in the Fig^c) and in Fig^a!), the E field is from upper left to the lower right coroner making an 
angle 30° with the horizontal, a semi-diagonal field (E field along 60° with the horizontal on the 
triangular lattice is not considered because it will lead to anisotropic current distribution). The 
presence of E field is shown by the long arrows. The B field is always perpendicular to the plane 
of the lattice and directed into the plane as indicated by the encircled dots. The descriptions given 
below are valid for all field configurations given in Fig|2 The black circles represent the occupied 
sites and the open circles represent the empty sites. The direction from which the central site is 
occupied is represented by a short thick arrow. In Fig^ the central site is always occupied from site 
2. The eligible empty sites for occupation due to E field are indicated by dotted arrows and thin 
solid arrows indicate the same due to B field. The number of eligible empty sites for occupation 
due to B field depends on the lattice structure. There are two such sites on the square lattice, as 
shown in Fig^a) and^fc), whereas there are three such sites available on the triangular lattice, 
as shown in Fig^c) and^d) by thin arrows. If the direction of B field is changed from 'into the 
plane' to 'out of the plane' only the sense of rotational constraint will be changed but the number 
of eligible empty sites for occupation will remain the same. However, the number of eligible empty 
sites for occupation due to E field depends on the direction of the E field. There is one such site 
for the horizontal E field on both the lattices whereas there are two such sites for the diagonal E 
field on the square lattice and for the semi-diagonal E field on the triangular lattice. Two eligible 
sites for occupation in case of diagonal or semi-diagonal E field orientation is the consequence of 
two equal components of the E field along those two nearest neighbours, as indicated by two dotted 
arrows m Figlflb) and FigflXd). 

After selecting the eligible empty sites for occupation, they are occupied with probability p. 
The coordinate of an occupied site in a cluster is denoted by (x,y). Periodic boundary conditions 
are applied in both directions and the coordinates of the occupied sites are adjusted accordingly 
whenever the boundary is crossed. At each time step the span of the cluster in the x and y directions 
L x = Xmax — x m i n and L y = y max — Vmin are determined. If either L x or L y of a cluster is found 
greater than or equal to L, the system size, the cluster is then considered as a spanning cluster. The 
critical percolation probability p c is the maximum probability below which there is no spanning 
cluster appears and at p = p c a spanning cluster appears for the first time in the system. 

3 Effect of Field Direction 

The critical properties of DSP model is studied numerically under diagonal E field on the square 
lattice and semi-diagonal E field on the triangular lattice. For the diagonal or semi-diagonal E 
field orientation, there are two components of the field on both the lattices along which the lattice 
sites could be occupied. Results obtained here are compared with that of already obtained results 
under horizontal E field orientations, single component of E field, on the respective lattices p3 Ej. 
Lattice size is varied from L = 2 7 to 2 11 in multiple of 2. Firstly the percolation thresholds (p c s) 
are determined for each field configuration on a given lattice. The cluster related quantities are 
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then evaluated at their respective p c s. Each data point is averaged over N to t = 5 x 10 4 clusters. 

3.1 Percolation threshold 

Percolation threshold p c is the maximum site occupation probability at which a spanning cluster 
appears for the first time in the system. p c s are determined here for the diagonal E field on the 
square lattice and semi-diagonal E field on the triangular lattice in the presence of a crossed B 
field. The probability to have a spanning cluster is given by P sp = n S p/Nt t where n sp is the number 
of spanning clusters out of Ntot clusters generated. The value of p c for a given system size L and a 
given E and B field configuration is determined from the maximum slope of the curve P sp versus 
p. Generally, p c depends on L for finite systems. The value of p c is then determined as function of 
L for a given lattice and field configuration. Data are then extrapolated to L — > oo. In FigEl p c {L) 
is plotted against the inverse system size 1/L for different field configurations on both the lattices. 
Circles represent p c s for the diagonal E field on the square lattice and diamonds represent the same 
for the semi-diagonal E field on the triangular lattice. The p c s of infinite systems are identified 
in the 1/L — » limit and marked by crosses. It is found that p c ~ 0.619 for diagonal E field on 
the square lattice and ~ 0.545 for semi-diagonal E field on the triangular lattice respectively. For 
the sake of comparison, p c s for the horizontal field configuration on both the lattices for different 
system size are also plotted in the same figure. The squares represent data for the horizontal E field 
on the square lattice and triangles represent the same for the horizontal E field on the triangular 
lattice. It was already obtained that p c ~ 0.655 1 and ~ 0.570.2] for the horizontal E field on 
the square and triangular lattices respectively as it is also shown here. As usual, the value of p c 
depends on the type of lattice as well as on the field situation. If the lattice structure is changed, 
the number of components of rotational constraint changes. If the direction of E field changes on a 
lattice, the number of components of directional constraint changes. The value of p c is decreased if 
the number of field components is increased either by the rotational constraint or by the directional 
constraint. It is expected, because, the clusters can grow rapidly if the number of field components 
are higher in a given lattice. 

3.2 Effective field and Spanning clusters 

Typical spanning clusters at p = p c for different field configurations on the square and triangular 
lattices of size L = 2 8 are shown in Fig|21 Spanning clusters generated on the square lattice for the 
horizontal and diagonal E field orientations are represented in Fig|2fa) and in Fig|3f6) respectively. 
FigEfc) and Fig|3{<i) represent the spanning clusters on the triangular lattice for the horizontal 
and semi-diagonal (30° with the horizontal) E field orientations respectively. The B field is always 
perpendicular to the E field and directed into the plane of the lattice as shown by encircled dots. 
The black dots indicate the external boundary of the cluster and the gray dots represent the 
interior of the cluster. It can be seen that the clusters are highly rarefied, anisotropic and have 
chiral dangling ends. Interestingly, the elongation of the clusters are not along the applied E field. 
This is due to the generation of the Hall field Eh perpendicular to the applied fields E and B. As 
a result, an effective field E e ff develops in the system and the elongation of the clusters are along 
the effective field. Since the direction of the effective field depends on the direction of the applied 
E field, the orientation of the cluster grown then depends on the direction of the applied E field. 
However, the clusters are not merely DP clusters along the effective field. Because, the presence of 
the rotational B field gives rise to chiral dangling ends at the boundary which has its own critical 
behaviour It is important to notice that the clusters are less anisotropic and more compact 
on the triangular lattice than on the square lattice. This is due to the extra flexibility given in 
the B field on the triangular lattice. However, the clusters look indifferent under different E field 
orientations on a given lattice except their direction of growth. 

3.3 Fractal dimension 

The fractal dimension df of the infinite or spanning clusters are determined by the box counting 
method. The number of boxes Nb(z) grows with the box size e as iYe(e) ~ e rf/ where df is the 
fractal dimension. In Fig0J Nb{c) is plotted against the box size e for the square (a) and triangular 
(6) lattices for all E field orientations. On the square lattice, the fractal dimensions are found as 
d f = 1.733 ± 0.005 for the horizontal E field and d f = 1.732 ± 0.006 for the diagonal E field. 
Similarly on the triangular lattice, it is found that df = 1.775 ±0.004 for the horizontal E field and 
df = 1.777 ± 0.005 for the semi-diagonal E field. The error is due to the least square fit. These 
estimates are also confirmed by finite size (FS) analysis. In FS analysis, df is determined employing 
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Soo ~ L d f , where is the mass of the spanning cluster. The FS estimates are df = 1.72 ±0.02 on 
the square lattice and df = 1.81 ±0.02 on the triangular lattice for the diagonal and semi-diagonal 
E field respectively. It is interesting to notice that the change in the direction of the E field has 
no effect on the values of df on both the lattices. However, the values of df are different on the 
square and triangular lattices as already reported in Ref.j^]. It was already found that df is higher 
on the triangular lattice than that of on the square lattice consistent with the observation that the 
spanning clusters are more compact on the triangular lattice than on the square lattice. For a given 
E field direction, the DSP model on these two lattices differ only by an extra rotational direction 
due to B field on the triangular lattice than that on the square lattice. The change in df on these 
two lattices is then only due to the different number of branching of the rotational constraint due 
to the B field in the presence of a directional field. The extra flexibility of the rotational constraint 
allows the clusters to penetrate more and more into itself on the triangular lattice than on the 
square lattice. As a result, the infinite clusters are less rarefied on the triangular lattice than on 
the square lattice. Since the clusters grow along the effective field direction, the orientation of 
the clusters is only changing with the orientation of the E field as the effective field direction is 
changing. Note that, the value of df obtained here is not only different from other percolation 
models but also smallest among that of OP (91/480), DP (1.765QII1) and SP (1.969HU). The 
spanning clusters are then most rarefied in DSP. Since df obtained in DSP is different from that of 
DP clusters, it could be inferred that the DSP clusters are just not DP clusters along the effective 
field. The DSP clusters then should have its own critical behaviour at the percolation threshold as 
it is already observed 1 . It is now important to check the effect of E field direction on the other 
critical exponents. 

3.4 Critical exponents and Scaling relations 

The values of some of the critical exponents are estimated and the scaling relations among them 
are verified in this section. A full scaling theory of DSP model is given in Ref.pQ. The cluster size 
distribution is defined as P s {p) = N s /Nt t where N s is the number of s-sited clusters in a total of 
Ntat clusters generated. Since the origin of a cluster is occupied with unit probability, the scaling 
function form of the cluster size distribution is assumed as 

P s (p) = s-^ 1 f[ S ^p- Pc )] (1) 

where t and a are two exponents. The form of P s {p) has already been verified and found appropriate 
for DSP ^H]. In order to verify the scaling relations among the critical exponents, the average 
cluster size, x = J2' s s Ps(p) ~ \p — J>c|~ 7 , and two other higher moments, \i = Z)s s2 Ps{p) ~ 
\p — p c \~ S and X2 = J2's s3 Ps(p) ~ \p — Pc\~ v , of the cluster size distribution P s (p) are measured 
generating finite clusters below p c . E field is applied diagonally on the square lattice and semi- 
diagonally on the triangular lattice in the presence of a crossed B field directed into the plane of 
the lattice. For a system size L = 2 , x is plotted against \p — p c \ for diagonal E field on the 
square lattice (circles) in FigEfa) and for semi-diagonal E field on the triangular lattice (diamond) 
in FigOJ&). Data corresponding to horizontal E field on the square (squares) and triangular 
(triangles) lattices are also plotted in the respective figures for comparison with the present 
data. It can be seen from Fig|S]that not only the absolute magnitude of the average cluster size x 
but also the scaling of \ with \p — p c | remains independent of the E field orientation on both the 
lattices. The value of 7 is found 1.87 ± 0.02 for the the diagonal E field whereas for the horizontal 
E field it was 1.85 ± 0.02 on the square lattice. Thus, the value of the exponent 7 is within error 
bar for both the E field orientations on the square lattice. A similar result is also obtained on 
the triangular lattice. It is found that 7 = 2.00 ± 0.02 for semi-diagonal E field whereas it was 
1.98 ± 0.02 for the horizontal E field on the triangular lattice. The errors quoted here include the 
least square fit error and the change in 7 with p c ± Ap c for Ap c = 0.0005. Similarly, estimates of 
the critical exponents 5 and 77 related to the higher moments of the cluster size distributions are 
also obtained. The observations are listed in Table ^ It can be seen that the values of all the 
critical exponents are very close and within error bars for different E field orientations on a given 
lattice. Finite size analysis of these critical exponents 7, 5 and In have also been made. The values 
of 7, 5 and r\ are determined on different system size L. In FigEj the exponents are plotted against 
the system size L. The extrapolated values of the exponents in the L — > 00 limit are marked by 
crosses. As expected, the finite size estimates of 7, 5 and r\ are converging to the Monte Carlo 
results on both the lattices as L — > 00. It is interesting to notice that the exponent values are 
close and within error bars for different E field orientation on a given lattice although they are 
significantly different (out of error bars) on the square and triangular lattices. The exponents 7, 5 
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and 77 satisfy a scaling relation 77 = 25 — 7. The scaling relation has already been verified in Ref.pQj 
and P] for horizontal E field and found holds true. The same scaling relation has also been tested 
with the obtained values of the exponents for diagonal/semi-diagonal orientations of E field and it 
is satisfied within error bars irrespective of E field orientations on both the lattices. 

Since the clusters generated here are anisotropic, there are then two connectivity lengths: £11 
along the elongation of the cluster and £j_ along the normal to the elongation. They are defined 
as q = 2E' s Rp)sPs(p)/E'sSPs(p) and £ = 2 £' s R\{s)sP s {p)/ £' s sP s (p) where fly (a) and 
R±(s) are the two radii of gyration with respect to two principal axes. The connectivity lengths, 
£|| ~ \P ~ Pc\~ u ^ and £j_ ~ \p — p c \~ u± diverge with two different critical exponents v\\ and uj_ as 
p —s- p c . The scaling behaviour of £y and £j_ are studied as function of p below p c for diagonal E 
field on the square lattice and for semi-diagonal E field on the triangular lattice. For lattice size 
L = 2 , £|| and are plotted against \p — p c \ in FigEta) for the square lattice and in FigEffr) for 
the triangular lattice for different E field orientations. The exponents z^y and v± are obtained as 

1.34 ± 0.02 and 1.13 =L 0.02 respectively for diagonal E field on the square lattice (circles). On the 
triangular lattice, they are found as v\\ = 1.37 ± 0.02 and v± = 1.24 ± 0.02 for the semi-diagonal 
E field (diamonds). Data for the horizontal E field are indicated by squares for the square lattice 
and by triangles for the triangular lattice. The values of U11 and u±_ are found with in error bars 
for different E field orientations on both the lattices, see Tabled The errors quoted here are the 
least square fit errors. Finite size estimates of these exponents are also made and it is found that 
they are consistent with that of the MC results. The exponent vn is almost same on the square and 
triangular lattices whereas u± on the triangular lattice is little higher than that of the square lattice 
value. Au = v\\—v±_ is ~ 0.21 for DSP on the square lattice whereas it is ~ 0.64 in DP^zj. The DSP 
clusters are then less anisotropic than that of DP clusters. This observation remains unchanged 
even if the direction of the effective field, along which the clusters grow, is changed. It is also 
noticed that the critical behaviour of the connectivity length £11 along the elongation of the cluster 
is similar to that of the connectivity length £ of ordinary percolation. The value of uu is almost 
equal to u{OP) = 4/3jH]. This is in good agreement with a recent study of magnetoresistance of a 
three-component composites of metallic film by Barabash et al^FQ. It is found that the hyperscaling 
relations 25 — 37 = (d— + u» and (d — dj)i , ± = 5 — 2^ are satisfied on the square lattice within 
error bars whereas they are satisfied marginally on the triangular lattice. It is already known that 
the hyperscaling is violated in DP|14j. 

The results obtained for different E field orientations on the square and triangular lattices are 
summarized in Tabled Comparing the data given in Table^h can be concluded that the values of 
the critical exponents and the scaling relations among them are independent of E field orientation 
or the number of components of E field irrespective of lattice structure. The direction of E field 
is only able to change the orientation of the clusters and not the critical properties. However, it 
has already been observed that the values of the critical exponents obtained in DSP model are 
not only different from that of the other percolation models but also are significantly different 
on the square and triangular lattices in 2D. It is already reported in Ref.j^] as a breakdown of 
universality between the square and triangular lattices. Microscopically, the DSP model on the 
square and triangular lattices differs only by an extra rotational direction due to B field on the 
triangular lattice than that on the square lattice for a given E orientation. The values of the 
critical exponents then depend on the number of branching due to the rotational field B but do 
not depend on the number of branching due to the directional field E. By reducing the number 
of branching of B field from three to two on the triangular lattice and keeping E field horizontal, 
an artificial model has been created and the critical exponents are re-estimated. The values of 
the critical exponents for this particular field combination are obtained as: df = 1.730 ± 0.005, 
7 = 1.87±0.01, 5 = 4.08±0.04 and rj = 6.34±0.08. Remarkably, the values of df and other critical 
exponents are found within error bars of that of the square lattice with horizontal E field and the 
un iver sality holds. It should be mentioned here that no breakdown of universality was observed in 
SP^J, percolation under only rotational constraint, between the square and triangular lattices. In 
DSP model, the presence of the directional constraint E on top of the rot atio nal constraint B may 
increase the number of spiral lattice trees. Since the spiral lattice trees|15] as well as the spiral 
self avoiding walks 16. have different scaling behaviour on the square and triangular lattices, the 
higher number of tree like structures in the clusters generated can induced a non-trivial effect on 
the critical properties of DSP at the percolation threshold and consequently the discrepancy in the 
exponent values on the two lattices could occur. 

Though there is a breakdown of universality in the cluster property of the DSP model on the 
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square and triangular lattices, it was already observed that the universalit y ho lds true for their hull 
(external perimeter) properties between the square and triangular lattices [T7]. The hull properties 
are also verified here by changing the applied E field direction and it is found that the hull properties 
are independent of the E field orientations. 

4 DSP Model under Variable Field Intensity 

So far, the DSP model has been investigated under constant external field intensities. It is also 
useful to understand the properties of disordered system as function of the field intensity. The 
intensities E and B of the external fields can vary here from to 1. It may be noted that for 
certain extreme values of the intensities E and B, the percolation model has four different varieties, 
DSP, DP, SP and OP. DSP corresponds to E = 1, B = 1; DP corresponds to E = 1, B = 0; SP 
corresponds to E = 0, B = 1 and OP corresponds to E = 0, B = 0. Thus, a complete phase 
diagram can also be obtained by studying the DSP model under variable field intensities between 
one and zero. If the field intensities are changed continuously from (1,1) to (0, 0) keeping E = B, 
it is expected to have a phase change from DSP to OP. Similarly, a phase change from DSP to SP 
can be obtained by changing E from 1 to keeping B = 1 and phase change from DSP to DP can 
be obtained by changing B from 1 to keeping E = 1. In order to study the effect of intensities of 
the external applied fields, one needs to generate percolation clusters under variable field intensities 
and characterize their geometrical properties at the percolation threshold. The intensity dependent 
percolation cluster could be generated by selecting nearest neighbours of an occupied site with 
appropriate probabilities, function of field intensities. To demonstrate the model, a square lattice 
of size L is considered. The E field is applied diagonally and the B field is applied perpendicular 
to the E field and directed into the plane of the lattice as shown in Fig^b). However, the results 
of this particular field configuration should be valid for all other field configurations. In order to 
occupy a site in the presence of external fields, a list of eligible sites for occupation is prepared first. 
The field intensities play role in the preparation of the list of eligible sites. If the field intensities 
are either E = 1 or B = 1 or both E = B = 1, the sites only in the favourable directions, defined 
by the directional and rotational constraints, are eligible for occupation. Sites in other directions 
are not considered for occupation in all three constrained percolation models DSP, DP and SP. The 
sites in the favourable directions of the fields are selected with probabilities proportional to the 
field intensities. If the field intensities are less than one, the sites in unfavourable directions of the 
fields will also be eligible for occupation due to scattering. A scattering field S is then introduced 
in the model. Sites in unfavourable directions are selected for occupation with probability S and 
it is defined as 

S = {I- E) x (1- B). (2) 

Notice that S = in all three cases: E = 1, B = 1 and E = B = 1. That means there is no 
scattering and consequently no sites in unfavourable directions are eligible for occupation as in the 
case of DP, SP or DSP. On the other hand, when E and B both are equal to zero then S = 1. 
There will be then no favourable direction. All empty nearest neighbours of an occupied site on 
the square lattice are selected with unit probability for occupation. The situation corresponds to 
OP. Thus, a change of phase from DSP to OP, DP or SP will be possible by changing E and B 
continuously. As soon as the list of eligible sites is prepared they are then occupied with percolation 
probability p. Below, the critical properties of DSP clusters are studied varying intensity of E and 
B fields. For each value of E and B, the percolation threshold p c is determined. Data are averaged 
over 10 5 clusters for each field intensity. 

5 Effect of field intensity 

The field intensities are changed in three different ways. First, E and B are changed from (1,1) to 
(0, 0) keeping E = B and accordingly DSP is expected to change over to OP. The corresponding 
data will be represented by circles. Second, keeping B = 1, E is changed from 1 to 0, then 
DSP should change to SP. These data are represented by squares. Third, keeping E = 1, B is 
changed from 1 to 0, thus DP will be reached from DSP. Triangles represent the corresponding 
data. Percolation threshold p c , at which a spanning cluster appears for the first time in the system, 
has been calculated as function of the field intensities E and B. Following the same procedure of 
section 3, p c s are identified as the percolation probability p corresponding to the maximum slope 
of the plot of spanning probability P sp against p. In Fig|Hl the values of p c s are plotted against E 
and B for a square lattice of size L = 1024. Each p c has been determined generating 10 5 spanning 
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clusters. It can be seen that as E and B goes to zero p c coincides with that of OP. Similarly, p c 
for DP and SP are also obtained for E = 1, B = and E = 0, B = 1 respectively. 

The percolation clusters are isotropic for OP and SP whereas they are anisotropic in the case of 
DP and DSP. For anisotropic clusters, there are two connectivity lengths £j_ and £ii and for isotropic 
clusters there is only one connectivity length. Since £j_ and £y diverge with two critical exponents 
u± and z/y respectively as p — > p c , the ratio of two connectivity lengths should diverge as 

\p — p c \~ Au where Au = vn — u±. For isotropic clusters, u\\ is equal to uj_ since the clusters grow 
equally in all directions. Thus, Au is zero for OP and SP clusters and Au is finite for anisotropic 
clusters of DP and DSP. In DSP, Au is approximately 0.21 pp| whereas in DP, it is approximately 
0.64 112] . The values of the connectivity exponents u± and v» are determined for different values of 
E and B. In Fig El Au is plotted against E and B. As the field intensities are taken to E = and 
B = or E = and B = 1 corresponding to OP and SP models the value of Au continuously goes 
to zero. As the field intensities are taken to E = 1 and B = corresponding to the DP model, 
the value of Au approaches to 0.64 as it is expected for DP. Not only the difference but also the 
absolute value of the connectivity exponent matches with the respective percolation models. Thus, 
there is a smooth crossover from DSP to OP, DP and SP at the appropriate values of E and B as 
they are changed continuously. 

In order to determine other critical properties, fractal dimension df of the spanning clusters and 
average cluster size exponent 7 are also determined, df is plotted in Fig^Jand 7 is plotted in FiglTTI 
against E and B. The value of df in different percolation model are: df{OP) = 91/48 ~ 1.896 9 , 
df(DP) « 1.765 Hm, d f (SP) « 1.957QH and d f (DSP) w 1.733 Q. The values of d f for OP, DP 
and SP are marked by crosses in FiglTUl It can be seen that the value of df changes continuously 
and approaches to the expected values of df of the respective percolation models as the intensities 
of E and B are changed appropriately. In Fig llll the values of 7 of differ ent percolation mo dels 
are also indicated by crosses, -/{OP) = 43/18 » 2.3890, -/(DP) » 2.2772QU1, -/(SP) « 2.19Q3] 
and -/(DSP) = 1.85 1 . Similar to the fractal dimension, the value of the critical exponent 7 is 
also changing continuously to that of the other percolation models as the intensities of E and B 
are changed appropriately. However, for most of the E and B values, the values of df and other 
critical exponents are found close to that of DSP. It could also be noticed that in the presence of 
small field intensities of E or B, the universality class of OP is changed. Similarly for DP and SP, 
a small intensity of the other field is able to change the universality class of the models to that of 
DSP. This is a new and important observation. 

A phase diagram for percolation models under external bias fields now can be plotted in the 
space of directional (E) and rotational (B) bia s fields. Different percolation models are identified 
on the E — B space by black circles in Fig ll2l Each point in this space corresponds to a second 
order phase transition point. The thick lines from left lower to right upper diagonal and parallel 
to the E and B axes of the phase diagram represent lines of second order phase transition points. 
The two ends of the diagonal line correspond to OP and DSP. The line parallel to E connects 
SP and DSP whereas the line parallel to B connects DP and DSP. The dotted lines define the 
regions of DP and SP as indicated. DSP model then can be considered as the most generalized 
model of percolation under external bias fields whose different limiting situations correspond to 
other percolation models like OP, DP or SP. 

6 Summary 

The directed spiral percolation model has been studied varying the applied field directions and 
their intensities on the square and triangular lattices in 2D. The critical properties of the cluster 
related quantities in this model were already found very different from the other percolation models 
like OP, DP and SP and accordingly DSP belongs to a new universality class. It is found here that 
the critical properties as well as the universality class of the DSP clusters remain invariant on the 
directions of the applied E and B fields. Change of E field direction corresponds to the change in 
the number of components of the E field on both the lattices. However, this is not the case for B 
field. Change in the number of components in the E field corresponds to the change in orientations 
of the clusters only whereas change in the number of components in the B field corresponds to 
the change in the universality class of DSP. Changing the intensity of the applied fields, a smooth 
crossover from DSP to OP, DP and SP is observed. A phase diagram is obtained for the percolation 
models under external bias fields. OP, DP or SP could be obtained as a limiting situation of DSP 
by changing the field intensities. The DSP model then can be considered as a most general model 
of percolation under external bias fields. The model will be applicable to the physical situations 
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where crossed directional and rotational force fields are present in disordered systems. 
Acknowledgment: SS thanks CSIR, India for financial support. 
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Square: 


TT '• ~i 1 I -J ■ 

Horizontal 1 : 


1.733 ± 0.005 


1.85 ± 0.01 


4.01 ± 0.04 


6.21 ± 0.08 


1.33 ± 0.01 


1.12 ± 0.03 




1.72 ±0.02 (FS) 














Diagonal: 


L.l 62 ± 0.0U6 


l.o f ± 0.01 


A r\/i _i_ n n/i 
4.U4 ± U.U4 


D.24 ± 0.08 


l.o4 ± U.U2 


i i q _t_ n no 
l.lo ± 0.02 






1.72 ±0.02 (FS) 












Triangular: 


Horizontal 0: 


1.775 ± 0.004 


1.98 ±0.02 


4.30 ± 0.04 


6.66 ±0.08 


1.36 ±0.02 


1.23 ±0.02 




1.80 ±0.03 (FS) 














Semi-diagonal: 


1.777 ±0.005 


2.00 ±0.02 


4.33 ± 0.04 


6.71 ±0.08 


1.37 ±0.02 


1.24 ±0.02 






1.81 ±0.02 (FS) 













Table 1: Numerical estimates of the critical exponents and the fractal dimension measured for the DSP clusters on the square and triangular 
lattices for different E field orientations. The values of the critical exponents and the fractal dimension are found independent of the orientation 
of E field on both the lattices. 
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Figure 1: Selection of empty nearest neighbors for occupation in a MC time step on the square 
(a, b) and triangular (c, d) lattices. Black circles are the occupied sites and the open circles are 
the empty sites. Thick long arrows represent directional constraint (E). The clockwise rotational 
constraint (B) is shown by encircled dots. The central site is occupied from site 2 and shown by 
short thick arrows. The eligible sites for occupation due to E field are shown by dotted arrows and 
thin solid arrows indicate the same due to B field on both the lattices. 
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Figure 2: Plot of critical probability p c (L) versus 1/L for different E field orientations. The 
symbols are: (□) for the horizontal E and (O) for the diagonal E field on the square lattice, (A) 
for the horizontal E and (O) for the semi-diagonal E field on the triangular lattice. The values of 
p c (L) are extrapolated to L — > oo and p c s for infinite systems are indicated by crosses. 
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Figure 3: Infinite spanning clusters at p = p c on the square (a, b) and triangular (c, d) lattices of 
size L = 2 8 . For the square lattice, the E field orientation is horizontal in (a) and it is diagonal 
in (b). For the triangular lattice, the E field orientation is horizontal in (c) and it is semi-diagonal 
(30° with the horizontal) in {d). For all E field configurations, the B field is directed into the plane 
of the lattice. The field lines are drawn from the origin of the clusters. Eh is the Hall field normal 
to both E and B. E e ff is the resultant field of E and Ejj- The gray dots represent the interior 
of the clusters and the black dots represent the cluster hull. The clusters are elongated along the 
effective field. 
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Figure 4: Number of boxes iVs(e) is plotted against the box size e for the square (a) and triangular 
(6) lattices to determine the fractal dimension df of the spanning clusters. Squares are for the 
horizontal E and circles are for the diagonal E field on the square lattice. Triangles are for the 
horizontal E and diamonds are for the semi-diagonal E field on the triangular lattice. The solid 
lines are guide to eye. The values of df are obtained as df = 1.732 ± 0.006 for diagonal E field on 
the square lattice and df = 1.777 ± 0.005 for the semi-diagonal E field on the triangular lattice. 
The value df is then independent of the E field orientation. 
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Figure 5: Plot of x against \p — p c \: (a) on the square and (6) on the triangular lattice for different 
E field configurations. The same symbol set of the previous figure is used. The solid lines are fitted 
to the data corresponding to the diagonal E field with slope 1.87 ±0.01 on the square lattice in (a) 
and to the semi-diagonal E field with slope 2.00 ± 0.02 on the triangular lattice in (6). 
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Figure 6: Plot of (a) 7, (b) 5 and (c) rj against 1/L on the square ( □, O) an d triangular (A, O) 
lattices for different E field configurations. The extrapolated values of the exponents are: 7 ~ 1.85 
and 1.87, 5 ~ 4.01 and 4.04 and ij ~ 6.21 and 6.24 for horizontal and diagonal E field on the square 
lattice. On the triangular lattice, the values are: 7 ~ 1.98 and 2.00, S rs 4.30 and 4.33 and r/ ~ 6.66 
and 6.71 for the horizontal and semi-diagonal E fields. The values of the critical exponents are 
within error bars for different E field orientations. 
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Figure 7: Plot of £11 and £± against \p — p c \: (a) on the square and (b) on the triangular lattice. 
The same symbol set of FigEJis used for different fields. The values of vn and v± are found with in 
error bars for different E field orientations. 
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Figure 8: Plot of percolation threshold p c against field intensities E and B for a square lattice of 
size L = 1024. For appropriate values of E and B, p c s correspond to that of respective percolation 
models. 




Figure 9: Plot of Av = z/y — v± against field intensities E and B. Av continuously goes to zero 
as the field intensities goes to E = B = for OP and B = 1, E = for SP respectively. Au is also 
approaching 0.64 corresponding to DP for E = 1 and B = 0. 




Figure 10: Plot of fractal dimension df of the spanning clusters against field intensities E and B. 
As E and B change continuously, the value of df approaches to that of the respective percolation 
models at the appropriate values of E and B. 
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Figure 11: Plot of average cluster size exponent 7 against field intensities E and B. 7 changes 
continuously to that of the respective percolation models at the appropriate values of E and B. 
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Figure 12: Phase diagram of the percolation models under crossed external bias fields. Different 
percolation models are represented by black circles. DSP changes continuously to OP, DP or SP 
for the field intensities of E and B are changed to the appropriate values. 
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